f=@(x,y,z)4*x.*z.*exp(-x.*x.*y-z.*z);
tic,
I=integral3(f,0,2,0,pi,0,pi,'RelTol',1e-20),
toc